Optimal purification of a generic n-qudit state 
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We propose a quantum algorithm for the purification of a generic mixed state p of a n-qudit 
system by using an ancillary n-qudit system. The algorithm is optimal in that (i) the number of 
ancillary qudits cannot be reduced, (ii) the number of parameters which determine the purification 
state |\E r ) exactly equals the number of degrees of freedom of p, and (iii) [<&) is easily determined 
from the density matrix p. Moreover, we introduce a quantum circuit in which the quantum gates 
are unitary transformations acting on a 2n-qudit system. These transformations are determined 
by parameters that can be tuned to generate, once the ancillary qudits are disregarded, any given 
mixed n-qudit state. 

PACS numbers: 03.67.-a, 03.67.Ac 
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I. INTRODUCTION 

Purification is one of the basic tools in quantum infor- 
mation science [1, 2]: Given a mixed quantum system S 
described by a density matrix p it is possible to intro- 
duce another ancillary system A, such that the state |\&) 
of the composite system S + A is a pure state and p is 
recovered after partial tracing over A: p = Tr^(|^')(4 f |). 
The ancillary system may be a physical environment that 
must be taken into account when doing experiments on 
the system S, but not necessarily so. It may be a fic- 
titious environment that allows us to prove interesting 
results about the system under investigation. 

Purification is a tool of great value in quantum infor- 
mation science, with countless applications, for instance 
in the study of the distance between quantum states [2] , 
of the geometry of quantum states [3] and of the quan- 
tum capacity of noisy quantum channels [4] . Besides its 
theoretical relevance, purification is interesting for exper- 
imental implementations of quantum information proto- 
cols requiring mixed state. While the direct generation 
of a mixture p of quantum states necessarily involves sta- 
tistical errors, this problem can be avoided if the purifi- 
cation state |^) is generated. Of course, the price to pay 
is that one has to work with the enlarged system S + A, 
including the ancillary system A. 

Due to the partial trace structure the purification state 
of a density matrix p cannot be uniquely denned, as 
any unitary transformation JJ a <8> Is acting non-trivially 
on the ancillary system only maps the state into a 
new state \^') = (Ua ® ls)|^) which is again a purifica- 
tion of p. In this paper, we propose a quantum protocol 
that selects a specific purification | vp) . Such purification 
turns out to be very convenient since the state \^) de- 



pends on a number of parameters exactly equal to the 
number of degrees of freedom of a generic mixed state p 
and is easily determined as a function of p. Furthermore, 
we design a quantum circuit given by a sequence of quan- 
tum gates acting on both the system and the ancillary 
qudits. The parameters which determine such quantum 
gates can be tuned to generate, once the ancillary qubits 
are disregarded, any mixed state p of system S. Finally, 
our protocol works for any system of n qudits and re- 
quires n ancillary qudits. We show that the number of 
ancillary qudits is optimal, that is, it cannot be reduced 
if we wish to design a quantum circuit capable of gener- 
ating any n-qudit mixed state. 

Our paper is organized as follows. To set the nota- 
tions, we first briefly review the concept of purification 
(Sec. II). Then we propose our purification scheme, mov- 
ing from simple to more and more complex cases. We 
start with the purification of a single qubit mixed state 
(Sec. Ill), then we proceed with the qutrit (Sec. IV) and 
the two-qubit (Sec. V) cases and finally we illustrate the 
generic n-qudit case (Sec. VI). Appendix A provides a 
short account of the diagrams of states, namely of a tool 
very useful for the purposes of the present paper. 



II. PURIFICATION 

Given a quantum system S described by the density 
matrix p, it is possible to introduce another ancillary 
system A, such that the state \^f) of the composite system 
is pure and 



p = Tr4| *><*!). 



(1) 



This procedure, known as purification, allows us to asso- 
ciate a pure state with a density matrix p. A generic 
pure state of the global system S + A is given by 
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with {|a)} and {\i}} basis sets for the Hilbert spaces Ha 
and Hs, of dimensions M and N, associated with the 
subsystems A and S. Given a generic density matrix for 
system S, 



N-l 
i,j=0 



(3) 



we say that the state |<3/) defined by Eq. (2) is a purifi- 
cation of p if 

M-l JV-1 

p = -& A (|*)<*|)= ^ ]T ^^-KXil- (4) 

a— i,j=0 

The equality between (3) and (4) implies 



Pij 



M-l 



r 1 r 1 * 

a=0 



(5) 



It is clear that (5) always admits a solution, provided 
the Hilbert space of system A is large enough. More 
precisely, it is sufficient to consider a system A whose 
Hilbert space dimension is the same as that of system S. 
Indeed, if we express the reduced density matrix using 
its diagonal representation, 



£pi|*><* 



a purification for the density matrix (6) is given by 



(6) 



(7) 



with orthonormal basis for Ha- This purification 

procedure requires the diagonalization of the density ma- 
trix p and therefore is in general only numerically feasible. 

In what follows, we propose a different purification 
scheme, which is optimal in that, for the purification 
of a generic n-qudit state, the number m = n of qu- 
dits of the ancillary system A cannot be reduced. While 
this is the case also for the well-known purification (7) 
based on the spectral decomposition (6), we anticipate 
that our method readily provides a purification state | <]/) 
that depends on a number of parameters exactly equal 
to the number of degrees of freedom of a generic p. Note 
that, even when the number of ancillary qudits is optimal 
(to = n so that M = N), the number 2N 2 - 2 of real 
coefficients C a i determining the purification state (2) is 
in general much larger than the number N 2 — 1 of real 
free parameters that must be set to determine a generic 
density matrix of size N. Different choices of the coef- 
ficients C a i are therefore possible. Our choice provides 
a purification state |\&) depending on a number N 2 — 1 
of real parameters exactly equal to the number of real 
freedoms of a generic mixed states p. Furthermore, the 
coefficient C a i in (2) can be easily determined from con- 
ditions (5), in spite of the fact that these equations are 



nonlinear. Finally, we will see in the next sections that 
our scheme suggest a very convenient quantum circuit for 
the preparation of a generic density matrix p. To illus- 
trate the working of our purification method, we discuss 
cases of increasing complexity, from a single qubit state 
to a generic n-qudit state. 



III. QUBIT 
A. Mixed-state purification 

We consider M = N = 2 and set 



Coi € M+, 

Cio € K+, Cn = 0. 



We first determine Cqi from (5): 



PiiHOnP + ICnPHCoil 2 , 



(8) 



(9) 



where the last equality follows from the second line of (8). 
Since we have also chosen Coi to be real and nonnegative, 
we simply obtain 

Coi = sfpTx- (10) 

Then we can determine Coo from the condition 

Poi = CooCqj + Ci C*i = CooCoi, (11) 

since Coi is already known from (9). Finally, knowing 
Coo, w e can derive Cio from the condition 

Poo = I Coo | 2 + |Cio| 2 . (12) 

Taking into account that Cio € M+, we have 

Cio - VPoo-ICool 2 . (13) 

Note that, in the special case in which p\\ = we can 
remove any ambiguity in the definition of the purification 
state I*) by setting Coo = 0. In this case, p = |0)(0| is 
already a pure state and its "purification" is l^) = |10). 
Alternatively, one can reshuffle the basis state according 
to |0) <-> 1 1) [8] . For a generic state p\\ ^ and the state 
|*) reads, in the {\ai) = |00), |01), |10), |11)} basis, as 
follows: 



Coo 
Coi 
Cio 
Cu 



Poi 



PoaPu — PwPoi 



Pn 




(14) 
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FIG. 1: Quantum circuit for the purification of a single qubit. 
The qubits run from top to bottom from the least significant 
(LSB) to the most significant (MSB). 



B. Mixed-state generation 
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FIG. 2: Diagram of states for the purification of a single qubit. 
Starting from the input state 1 00) , information flows on the 
thick lines. We esplicitly indicate at the right hand side that 
the coefficientes Coi and Cio are real, while Cn = 0. 



In this subsection, we provide a quantum circuit gen- 
crating the state \^f) of (14), namely the purification of a 
generic single-qubit density matrix p. After disregarding 
the ancillary qubit (this corresponds to performing the 
partial trace of the density matrix |^')( V E'| over the ancil- 
lary qubit), we obtain the mixed state p. Therefore, we 
end up with an experimentally viable procedure for the 
generation of a generic mixed single-qubit state by means 
of a two-qubit state subjected??? to controlled unitary 
transformations. 

The quantum circuit generating the state is shown 
in Fig. 1. A square box with a greek letter inside (here, 
a or 9) stands for a rotation operator. Its matrix repre- 
sentation in the {|0), |1)} reads as follows: 



R(a) 



cos a — sm a 
sin a cos a 



(15) 



The full circle with —<j> above is the phase-shift gate, 
defined by the diagonal matrix 



PHASE(- 



diag(l,e 



(16) 



As an overall phase factor is arbitrary, the action 
of this gate is equivalently represented by the matrix 
diag(e l< ^, 1). In the controlled-gates, the empty circle on 
the control qubit means that the gate acts non trivially 
(differently from identity) on the target qubit if and only 
if the state of the control qubit is |0). 

On the other hand, any single-qubit density matrix can 
be generated by means of the quantum circuit in Fig. 1. 
Given the input state = |00), the output state is 



1*/) = 



cos a cos ve 
cos a sin ( 
sin a 




(17) 



This state is equal to the purification state (14), provided 
we set 



sin a = Cio, 
cos a sin 9 = Cn , 
cos a cos = Cot- 



(18) 



From the first equation we determine a, then from the 
second 9 and finally from the third (j), as a function of 
the coefficients C a i , which in turn are determined by our 
purification protocol from the density matrix elements 
Pij. Note that, since by construction Coi, Cio > 0, we 
can take a, 9 £ [0, ?] . Finally, the phase <f> £ [0, 2tt). 

For a straightforward extension of the results presented 
in this section from the purification of a single qubit 
to more complex systems it is convenient to express 
the quantum circuit in Fig. 1 in terms of diagrams of 
states [5], of which a very brief account is given in Ap- 
pendix A. The diagram of states corresponding to the 
purification circuit of Fig. 1 is shown in Fig. 2. Note 
that, given the input state |00), the output state (17) is 
immediately written following the information flow along 
the thick lines of the diagram of states. We stress that 
other optimal purifications, where the purification state 
is determined by N 2 — 1 = 3 real parameters, are possi- 
ble. Our choice corresponds to setting Cn = and Coi, 
Cio rea l- We will see that the diagrams of states imme- 
diately lead to optimal purifications also for arbitrarily 
complex systems. 



C. "Invasion" of the Bloch ball 

Using the quantum circuit in Fig. 1 we can write a 
generic single-qubit state as 



P = z2PkHk){i>k\ 
fe=i 



(19) 



with pi — cos 2 a, p2 = sin 2 a, \ipi) generic single-qubit 
pure state and \ip2) = |0). The matrix representation of 
p in the {|0), |1)} basis is given by 



p = cos a 
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cos#sin0e~ 
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(20) 
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FIG. 3: "Invasion" of the Bloch ball: cos 2 a = (top left), fjj 
(top right), 2^ (middle left), 3^ (middle right), 4^ (bottom 
left), f (bottom right). 



IV. QUTRIT 



We consider M = N = 3 and set 



C02 G R+, 

Cii e R+, C12 = 0, 

C20 € K+, C21 = C22 



0. 



(21) 



We first obtain C02 from 

P22 = \c 02 \ 2 + \c 12 \ 2 



I C22 1 2 — I C02 1 2 



Once C02 > is determined as Ci 



02 



(22) 
P22, we obtain 



C01 and Coo from 

P12 = C01C02 + CiiC7 12 + C21C22 = C01C02, (23) 

P02 = CooCo 2 + CloC 12 + C2oC2 2 = CooCo2- (24) 

Then we obtain C\\ from 

Pxx = |Coi| 2 + |C u | 2 + |C 21 | 2 = |C 01 | 2 + C 1 2 1 , (25) 
and, finally, C10 and C20 from 

POl = CooCq! + CiqC^ + C20C21 = CooC*^ + CioCn, 



Poo 



IC00I 2 



ICxnl 2 



20- 



(26) 
(27) 



We stress that conditions (21) lead to a purification 
state determined by iV 2 — 1 = 8 free real parameters, 
exactly corresponding to the number of real freedoms 
needed to set a generic density matrix for a qutrit. Con- 
ditions (21) are readily derived if the purification of a 
generic qutrit state is implemented by means of the di- 
agram of states shown in Fig. 4. In this figure, a box 
with two geek letters written on top of it (for instance, 
oti and a 2 ) represents a unitary transformation whose 
matrix representation has, in the {|0), |1), |2)} basis, the 
first column given by 



cos OL\ COS OL2 

cos a.\ sin c*2 
sin a.\ 



(28) 



The last equality corresponds to the usual Bloch-ball rep- 
resentation of the (generally mixed) single-qubit states. 
It is clear that, once a is fixed, Eq. (20) represents a 
surface in the (X, Y, Z)-space, obtained after contracting 
the pure-states (unit radius) Bloch-sphere of the factor 
cos 2 a and translating it in the positive Z-direction by 
sin 2 a. Plots of this surface for different values of a are 
shown in Fig. 3. It is clear that all the points of the Bloch 
ball are recovered when a goes from to ^. We can say 
that there is an "invasion" of the Bloch ball starting from 
the north pole. 



Such transformation maps the input state |0) into 

cosai cosa2|0) + cosai sina 2 |l) + sinai|2). (29) 

Finally, the box with the letter 6>3 on top of it represents 
the rotation R(03) acting on the two-dimensional sub- 
space spanned by the state |0) and |1), with R defined 
by (15). 

Given the input state |00), the output purification state 
|^) = ; C a i\ai) can be immediately written by fol- 
lowing the thick lines of the diagram of states in Fig. 4. 
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FIG. 4: Diagram of states for the purification of a single 
qutrit. To simplify the plot, only the thick lines correspond- 
ing to the information flow are shown inside the boxes. The 
angles a t , 9j G [0, , while the phases <j>k G [0, 2tt). 



We obtain 



f Coo 


= COS OL\ COS Of 2 COS 


9\ cos 


9 2 e 




= COS OL\ COS Of 2 COS 


9\ sin ( 


he' 




= cos a.\ cos a 2 sin 


?i, 




Cio 


= cosai sina2 cos 


9 3 e* 




Cn 


= cos a\ sin a2 sin ( 






Cl2 


= o, 






C20 


= sin a.\ , 






C2I 


= 0, 






< C22 


= 0. 







These relations can be easily inverted to obtain the pa- 
rameters {aj, 8k,4>i} as a function of the coefficients C a % 
and, therefore, of the elements of the density matrix p. 
Therefore, the quantum circuit represented by the dia- 
gram of states of Fig. 4 can be used to generate any given 
qutrit state p, once the ancillary qutrit is disregarded. Fi- 
nally, we point out that the number of real parameters 
{cti,9k,4>i} that determine the state \^f) is equal to 8, 
that is, exactly to the number of parameters needed to 
determine a (generally mixed) single-qutrit state p. 

It is clear from the purification drawn in Fig. 4 that 
we can write a generic single-qutrit state as 

3 

P = ^2pk\ipk)(ipk\, (31) 

k=l 

with pi — cos 2 oti cos 2 a2, P2 — cos 2 ai sin 2 CV2, P3 — 
sin ai, l^i) generic single-qutrit pure state, 1-02) pure 
state residing in the two-dimensional subspace spanned 
by |0) and |1), and \tjj3) = |0). The picture developed in 
Sec. Ill C about the "invasion" of the single-qubit Bloch 



ball by means of suitably scaled and translated pure- 
qubit Bloch spheres may be generalized to the single- 
qutrit case. The role of the Bloch sphere is here played 
by the surface of single-qutrit pure states and the volume 
of all single-qutrit states is "invaded" when the parame- 
ters pk are varied, with the constraint ^2 k Pk = 1< 

V. TWO QUBITS 

We consider M = N = 4 and set 
C03 € K+, 

^ C\2 G K+, C13 = 0, ,^2% 

C21 G K.+, C22 = C23 = 0, 
k C30 G R+, C31 = C32 = C33 = 0. 

We first obtain C03 from 

P33 = |C 03 | 2 + ICial 2 + I C 23 1 2 + |C 33 | 2 = I C 03 1 2 , (33) 
then C02, Cox-, and Coo from 

P23 = C02C03 + Ci2C* 3 + C22CJ3 + C32CI3 = C02C03, 

(34) 

Pl3 = CoiCq 3 + CiiC* 3 + C21C23 + C31CI3 = C01C03, 

(35) 

P03 = CooCg 3 + CioC* 3 + C20C23 + C30C33 = C00C03, 

(36) 

then C12 from 

P22 = |Co2| 2 +|C 12 | 2 +|C22| 2 +|C32| 2 = |C 02 | 2 +C 2 2 , (37) 
then C\\ and C\o from 

P\2 = CoiCq 2 + C11C12 + C2lC 22 + C31C32 

= CoiCq 2 + CiiCi 2 , (38) 

P02 = CooCq 2 + C10CJ2 + C20C22 + C30C32 

= C00C02 + C\oC\2- (39) 

then C21 from 

P11 = |C i| 2 + |C 11 | 2 + |C 2 i| 2 + |C 3 i| 2 

= |C i| 2 + |C 11 | 2 + C 2 2 1 , (40) 

then C20 from 

P01 = CooCqi + C\qC 1x + C20C21 + C^oC^ 

= CooCqi + CioC^ + C20C21, (41) 

and, finally, C30 from 

Poo = |Coo| 2 + Idol 2 + |C 20 | 2 + Clo. (42) 

As for the previous examples, we point out that the 
number ./V 2 — 1 — 15 of free parameters determining the 
purification state cannot be reduced given a generic 
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FIG. 5: Diagram of states for the purification of a two-qubit 
state. To simplify the plot, only the thick lines corresponding 
to the information flow are shown inside the boxes. The angles 
cti, 6j € [0, , while the phases cifc E [0, 2tv). 



two-qubit mixed state p. Moreover, conditions (32) are 
determined from the diagram of states for the purifica- 
tion of a generic two-qubit state shown in Fig. 5. Fol- 
lowing the information flow from the input states 1 0000) 
we can immediately write down the output purification 
state \SS?) = ^2 a i C a i\ai). We obtain 



Coooo 




COS OL\ COS OL2 cos 


9\ cos 


9ze l 


Coooi 




COS OL\ COS «2 cos 


9 1 sin 


he 1 


Cboio 




cos ot\ cos a2 sin 


9i cos 


9 4 e l 


Coon 




cos tt\ cos (22 sin 


9i sin t 


h, 


Coioo 




cos ai sin oli cos 


?2 cos 


%e l 


Choi 




cos ot\ sin a,i cos 


9 2 sin t 


he 11 


Cono 




cos ot\ sin a-i sin t 


h, 




< Com 




o, 


9 6 e* 




Ciooo 




sin a.\ cos c*3 cos 




Ciooi 




sin a.\ cos a3 sin t 






Cioio 




o, 






Cion 




o. 






Cnoo 




sinai sincv3, 






Cnoi 




o, 






Cino 




0. 






k Cini 




0. 







As in the previous cases, we can invert these equations 
and determine the parameters {oej, 9k,4*i} m terms of the 
coefficients C a i and, therefore, of the elements of the den- 
sity matrix p. 

We can see from Fig. 5 that a generic two-qubit state 



can be written as 

4 

P = 5>*NfcXM (44) 
fc=i 

with pi — cos 2 ol\ cos 2 cx.2, Vi — cos 2 ai sin 2 a2, P3 = 
sin 2 a.\ cos 2 Ck3, and p4 = sin 2 ot\ sin 2 a?,. The "invasion" 
picture developed in Sec. IIIC can be extended also to 
the present case, with the volume of all two-qubit states 
"invaded" when the parameters pk are varied, under the 
constraint ^2 k Pk = 1- Note that the number of real pa- 
rameters {ai,^,^} used to determine the state |\&) is 
15, that is, exactly the number of parameters required to 
determine a generic two-qubit state. 

VI. n QUDITS 

We consider M = N = d n and set 

Co,AT-l S R+, 

Ci,jv-2 S R+, Ci^jv-i = 

C2.N-3 € K+, C2,Ar-2 = C2.N-I — 0, 

< . ' (45) 

Cat-i^o S R+, CiV-1,1 = Cjy-1,2 = • • • 
= Cpf-1,N-1 = 0. 

The general procedure for determining the coefficient 
C a i is clear from the previous examples. 

• We first determine Co,jv-i from pn-i,n-i, 

• then Coj from /Oj,jv— l; with j = N — 2, 0, 

• then Ci i AT_ 2 from PN-2.N-2, 

• then C\j from Pj,N-2, with j = N — 3, 0, 

• and, finally, Cv-i^o from poo- 

This purification is optimal as the number of qudits 
of the ancillary system cannot be reduced. The number 
of real free parameters that must be set to determine a 
density matrix of size d n is d 2n — 1 [the —1 term comes 
from the normalization condition Tr(p) = 1]. There- 
fore, an ancillary system of n — 1 qudits is not sufficient 
to purify p, as a pure state |^) in the Hilbert space of 
n + (n — 1) = 2n — 1 qudits has 2d 2 ™ -1 — 2 freedoms (the 
—2 term is due to the normalization condition for \^} and 
to the fact that a global phase factor in |\&) is arbitrary). 
This number is not sufficient as 2d 2n ~ 1 — 2 < d 2n — 1 
for any d > 2, n > 1. On the other hand, as illustrated 
Fig. 6, the number of free real parameters in our purifi- 
cation method is exactly d 2n — 1. From the schematic 
drawing in Fig. 5 we can also see that the n-qudit state 
can be written as 

d n 

P = X>fch^>vM (46) 
fe=i 
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is optimal, in that the number n of ancillary qudits used 
for the purification cannot be further reduced. Moreover, 
our algorithm can also be seen as a quantum protocol for 
the generation of a generic n-qudit state by means of suit- 
able unitary operations applied both to the system and to 
the ancillary qudits, with the ancillary qudits eventually 
disregarded. While also the well-known purification (7) 
based on the spectral decomposition (6) uses n ancillary 
qudits, our method is optimal in that it readily provides 
a purification state that depends on a number of param- 
eters exactly equal to the number of degrees of freedom 
of the generic n-qubit state that we wish to purify. 

It is well known that the purification |\&) of a generic 
mixed state p cannot be uniquely determined, as the par- 
tial trace over the ancillary qudits is invariant under any 
unitary transformation Ua <S> Is acting non trivially on 
the ancillary qudits only. Indeed, we have 



p = TMI*}<*|) = TMI*'X*'I 



(48) 



FIG. 6: Schetch of the diagram of states for the purification of 
a n-qudit state by means of 2n qudits. The diagram has d 2n 
lines, one for each state of the computational basis. We can 
cluster groups of d n lines. The constraints and the number of 
real freedoms for each cluster are highlighted. If we add these 
numbers plus the d n — 1 independent weights in the mixture 
(46), we obtain a total number of d 2n — 1 degrees of freedom. 



with the pure states from \ipi) to \ipd n ) residing in sub- 
spaces of decreasing dimension, from d n to 1 . 

We note that the tensor product structure of many- 
qudit quantum systems does not play any role in our 
purification protocol. That is to say, the same purifica- 
tion scheme applies to a system of n qudits or to a single 
system of size N = d n . Of course, the practical im- 
plementation of the protocol will depend on the specific 
quantum hardware at disposal. 

From the mathematical viewpoint, our purification 
protocol can be seen as the Cholesky decomposition [6] 
of the density matrix p. If we consider the coefficients 
C a i as elements of a d n x d n matrix C, then its transpose 
D = C T is a upper triangular matrix and Eq. (5) reads 

p = C T C* = DD T * = DD\ (47) 

namely it is the Cholesky decomposition of the density 
matrix p. Such decomposition is unique when p is pos- 
itive, while the ambiguities arising when p is singular 
may be removed by suitable prescriptions or reshuffling 
as previously discussed for the single-qubit case. 

It is interesting to remark that the Cholesky decompo- 
sition has already been used in the context of quantum 
information science in parametrizing the density opera- 
tor in order to guarantee positivity [7]. The purpose of 
Ref. [7] was to present a universal technique for quantum- 
state estimation. 

VII. FINAL REMARKS 

In summary, we have proposed an algorithm for the 
purification of a generic n-qudit state. This algorithm 



with 



|*'> = (U A ®t s )\^) 



(49) 



It is sufficient to add the unitary transformation Ua <8> Is 
at the end of our purification protocol to obtain any 2n- 
qudit purification (SP') of a n-qudit state p. On the other 
hand, we can say that our protocol selects a very conve- 
nient purification as the coefficients of the wave function 
|\E'} in the computational basis are easily determined from 
the density matrix p and the quantum circuit generating 
p can be immediately drawn. 



APPENDIX A: DIAGRAMS OF STATES 

Diagrams of states [5] graphically represent how quan- 
tum information is elaborated during the execution of a 
quantum circuit. In the usual way of drawing a quan- 
tum circuit [1, 2] each horizontal line represents a qubit. 
In contrast, in diagrams of states we draw a horizontal 
line for each state of the computational basis. Therefore, 
diagrams of states are less synthetic but may help us to 
clearly visualize quantum information flow in a quantum 
circuit. 

For the purposes of the present paper, it will be suffi- 
cient to show the diagrams of states for elementary single- 
qubit quantum gates. The phase-shift gate PHASE(0), 
defined by Eq. (16), and the rotation gate R(0), defined 
by Eq. (15), are shown in Fig. 7 and Fig. 8, respec- 
tively. Finally, the generation of a generic single-qubit 
state cos#|0) + sin6*e l *|l) starting from the input state 
|0) is shown in Fig. 9. The information flows on the thick 
lines, from left to right, while thinner lines correspond to 
absence of information. Note that, following the thick 
lines, the final state can here be immediately written. 
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FIG. 7: Quantum circuit (left) and diagram of states (right) 
for the phase-shift gate. 




FIG. 8: Quantum circuit (left) and diagram of states (right) 
for the R(6) gate. 




FIG. 9: Quantum circuit (left) and diagram of states (right) 
for the generation of a generic single-qubit state. To simplify 
the plot, only the thick lines corresponding to the information 
flow are shown inside the boxes. 
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